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Abstract 



We present a measurement of the mass difference Anid for the meson and 
the statistical power of the b flavor tagging methods used. The measurement 
uses 90 pb~^ of data from pp collisions at ^/s = 1.8 TeV collected with the 
CDF detector. An inclusive lepton trigger is used to collect a large sample 
of B hadron semileptonic decays. The mass difference Am^ is determined 
from the proper time dependence of the fraction of B hadrons that undergo 
flavor oscillations. The flavor at decay is inferred from the charge of the 
lepton from semileptonic B decay. The initial flavor is inferred by determining 
the flavor of the other B hadron produced in the collision, either from its 
semileptonic decay (soft-lepton tag) or from its jet charge. The measurement 
yields Am^ = 0.500 ± 0.052 ± 0.043 ^ps~^, where the first uncertainty is 
statistical and the second uncertainty is systematic. The statistical powers 
(eZ)^) of the soft-lepton and jet-charge flavor taggers are (0.91±0.10±0.11)% 
and (0.78 ± 0.12 ± 0.08)%, respectively 

PACS numbers: 12.15.Ff, 13.20.He, 14.40.Nd 



4 



Typeset using REVT^ 



I. INTRODUCTION 



In the Standard Model of electroweak interactions ||T]|, the quark mass eigenstates are 
related to the weak eigenstates via the unitary 3x3 Cabibbo-Kobayashi-Maskawa (CKM) 
Matrix Vckm 0- The nine elements in this matrix, Vij, where i = u, c, t and j = d, s, b, are 
completely determined from three angles and a phase. A nonzero phase gives CP violation in 
the weak interaction. Measurements of decays of hadrons containing b quarks {B hadrons) 
are of great interest because they determine the magnitudes of five of the nine elements of 
Vckm as well as the phase. 



A. The Unitarity Triangle 

The unitarity of Vckm leads to nine unitarity relationships, one of which is of particular 
interest : 

v^dV;; + KdK; + Vt,v,i = 0. (1.1) 

In the complex-plane, this sum of three complex numbers is a triangle, commonly referred to 
as the Unitarity Triangle. Measurements of the weak decays of B hadrons and the already 
known CKM matrix elements determine the magnitudes of the three sides of the Unitarity 
Triangle, and CP asymmetries in B meson decays determine the three angles. The primary 
goal of B physics in the next decade is to measure precisely both the sides and angles of 
this triangle and test consistency within the Standard Model. 



We can use several approximations to express eq. |1.1| in a more convenient form. The 
elements V^d — 1 and V^d — —A = — sin 9c, where 9c is the Cabibbo angle, are well measured. 
Although the elements Vtb and Vts are not well measured, the theoretical expectations are 



that Vtb — I and Vts — —Vet- With these assumptions, eq. |LT| becomes 



V* V 

_i_J^ = o. (1.2) 



Measurement of Am^, the subject of this paper, directly impacts the determination of 
Vtd/Vs. 



B. Determining Vtd and Vts From Neutral B Meson Flavor Oscillations 

Second-order weak processes transform a neutral B meson into its antiparticle: B^ ^ B^ ^ 
giving a probability for a B^ to decay as a B^ that oscillates with time. The frequency of 
these oscillations is the mass difference Am^ between the B mass eigenstates, which are 
linear combinations of the flavor eigenstates B^ and B^. The mass difference is proportional 
to iVJfoVfdP, so in principle, a measurement of Am^ determines this product of CKM ma- 
trix elements. In practice, however, large theoretical uncertainties limit the precision of 
Vtd- The same problem exists in determining Vts from Am^, the frequency of 5^ — >■ 5° 
oscillations. These theoretical uncertainties are reduced in determining the ratio |VJd/Vtsp 
from Arrid/ Arris- Unfortunately, at this time, attempts to measure Arris have only led to 
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lower limits. The determination of Am^ is a key future measurement of B hadrons, since 
l^d/Visp determines the magnitude of one of the sides of the Unitarity Triangle, as expressed 
in eq. |L^. 

Measurements of Am^ and Am^ require determining the initial flavor of the B meson, 
that is, whether the B meson contained a b quark or a 6 antiquark. Flavor determination 
is also crucial in the measurement of CP violation in the decays of neutral B mesons to 
CP eigenstates. 

This paper describes a measurement of Am^ using data collected by the CDF experi- 
ment from pp collisions with a center-of-mass energy of 1.8 TeV produced by the Fermilab 
Tevatron. In addition, the measurement of Am^ is used to demonstrate the performance of 
two methods (described below) of identifying the flavor of B hadrons in the environment of 
pp collisions. These methods of flavor identification will be important in the measurement 
of CP violation in the decays of neutral B mesons and in the study of Am^ 0. 

C. Previous Measurements of Am^ from CDF 

The CDF collaboration has exploited the large b quark cross-section at the Fermilab 
Tevatron to make several precision measurements of the properties of B hadrons, including 
lifetimes and masses [^]. Although the b quark production cross-section at the Tevatron 
is large, the pp inelastic cross-section is three orders of magnitude larger, so specialized 
triggers are required to collect large samples of B hadrons. To date, the triggers that have 
been utilized are based on leptons (electrons and muons) or dileptons. Some analyses use the 
semileptonic decays of B hadrons, B — > iuX; some use the semileptonic decays of charmed 
particles from B hadron decay {e.g. B DX, followed hj D ^ iuY); and some use leptonic 
J/ip decays {B J/ipX, followed hj J/ip ^ 

Previous measurements of Am^ from CDF were based on data samples collected with 
a dimuon trigger p and single-lepton triggers (£ = e,/i) The single-lepton triggers 
were used to partially reconstruct approximately 6,000 B^ mesons via their semileptonic 
decays 5° ^ e+uD^*^-X (in this paper, reference to a particular decay sequence implies 
the charge-conjugate sequence as well). The analysis reported in this paper uses the same 
data sample collected with this trigger, but increases the number of B mesons by over an 
order of magnitude by inclusively reconstructing B hadrons that decay semileptonically. The 
inclusive reconstruction is made possible by the relatively long lifetime of B hadrons: the 
decay point of the B hadron is typically separated from the production point (the primary 
vertex) by a couple of millimeters. The inclusive reconstruction is based on identifying 
this decay point by associating the trigger lepton with other charged decay products to 
reconstruct a secondary vertex. 

D. Method of Measuring Am^ 

The oscillation frequency Am^ can be found from either a time independent measurement 
(that is, from the total number of B^^s that decay as 5°'s) or from a time dependent 
measurement (that is, from the rates that a state that is pure B^ at t = decays as either a 
B^ or 5" as a function of proper decay time t). The latter technique has better sensitivity 
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and allows a simultaneous study of the tagging methods since the amplitude of the oscillation 
depends on the effectiveness of the tagging method. 
The expected rate is 

r{B° S°) = — e-*/"^(l + cos(Amdt)); (1.3) 

V{B° B^) = — e-*/^^(l - cos(Amdt)), (1.4) 

where tb is the mean lifetime of the two mass eigenstates of the B^ (the difference in lifetime 
of these two eigenstates is very small and has been neglected), and t is the proper decay 
time of the B^ in its restframe. To measure this rate, we need to make three measurements: 
(1) the proper decay time, (2) the B flavor at decay, and (3) the produced B flavor. We 
determine the proper time by measuring the distance from the production point to the decay 
point in the laboratory frame combined with an estimate of the 5° momentum. The flavor 
at decay is determined from the charge of the lepton, assuming it comes from semileptonic 
B decay. In our measurement of Am^ using B^ i^uD^*^" X decays, the flavor at 
production was identified using a same-side tagging technique based on the electric charge 
of particles produced in association with the B^. This flavor tag has also been applied ^ to 
a sample of B^/B^ —>■ J/ipKg decays to measure the CP asymmetry. To identify the flavor 
at production in this analysis, we rely on the fact that the dominant production mechanisms 
of b quarks in pp collisions produce bb pairs. The flavors of the B hadrons are assumed to 
be opposite at the time of production. In this paper, we identify the flavor of the other 
B hadron using two techniques: the soft-lepton flavor tag (SLT) and the jet-charge flavor 
tag (JCT). 

Several precise measurements of Am^ have been published by experiments operating 
on the T{4S) and resonances. The measurement of Am^ presented here is competitive 
in precision and in addition quantifies the performance of the flavor tags, which are crucial 
for future measurements of CP violation in the decays of B mesons and the measurement 
of Arris at a hadron collider. The jet-charge flavor tag is a powerful technique in studies of 
neutral B meson flavor oscillations by experiments operating on the Z° resonance ||10[ . This 



analysis is the first application of the jet- charge flavor tag in the environment of a hadron 
collider. 



II. THE CDF DETECTOR 

The data sample used in this analysis was collected from 90 pb^^ of pp collisions recorded 
with the CDF detector at the Fermilab Tevatron. The CDF detector is described in detail 
elsewhere |ll|Jl^ . We summarize here the features of the detector that are important for this 



analysis. The CDF coordinate system has the z axis pointing along the proton momentum, 
with the X axis located in the horizontal plane of the Tevatron storage ring, pointing radially 
outward, so that the y axis points up. The coordinates r, 0, and 9 are the standard cylindrical 
coordinates. 

The CDF spectrometer consists of three separate detectors for tracking charged particles: 
the silicon vertex detector (SVX), the vertex detector (VTX), and the central tracking 
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chamber (CTC), which are immersed in a magnetic field of 1.4 Tesla pointed along the 



+z axis. The SVX ||T3[ consists of four concentric cylinders of single-sided silicon strip 
detectors positioned at radii between 3 cm and 8 cm from the beam line. The strips are 
oriented parallel to the beam axis and have a pitch of 60 /im in the inner three layers 
and 55 /xm on the outermost layer. The SVX is surrounded by the VTX, which is used 
to determine the z coordinate of the pp interaction (the primary vertex). Surrounding 
the SVX and VTX is the CTC. The CTC is a drift chamber that is 3.2 m long with 84 
layers of sense wires located between a radius of 31 cm and 133 cm. The sense wires are 
organized into five axial superlayers and four stereo superlayers with a stereo angle of 3°. The 
momentum resolution of the spectrometer is Spx/pr = [(0.0009GeV/c ■ px)"^ + (0.0066)^] 2, 
where pt is the component of momentum transverse to the z axis {pt = p ■ sin^). Charged 
particle trajectories reconstructed in the CTC that are matched to strip-clusters in the SVX 
have an impact parameter resolution of 6do = [13 -|- (40 GeV/c)/^^] /^m, where the impact 
parameter do is the distance of closest approach of the trajectory to the beam axis in the 
plane perpendicular to the beam axis. The outer 54 layers of the CTC are instrumented to 
record the ionization (dE/dx) of charged tracks. 

Surrounding the CTC are the central electromagnetic calorimeter (CEM) and the cen- 
tral hadronic calorimeter (CHA). The CEM has strip chambers (CES) positioned at shower 
maximum, and a preshower detector (CPR) located at a depth of one radiation length. Be- 
yond the central calorimeters lie two sets of muon detectors. To reach these two detectors, 
particles produced at the primary vertex with a polar angle of 90° must traverse material to- 
taling 5.4 and 8.4 pion interaction lengths, respectively. The trigger system consists of three 
levels: the first two levels are implemented in hardware. The third level consists of soft- 
ware reconstruction algorithms that reconstruct the data, including three-dimensional track 
reconstruction in the CTC using a fast algorithm that is efficient only for p^ > 1.4 GeV/c. 



III. DATA SAMPLE SELECTION 

The sample selection begins with data from the inclusive e and fi triggers. At Level 2, 
both of these triggers require a track with pt > 7.5 GeV/c found by the central fast tracker 
(CFT) [|14|, a hardware track processor that uses fast timing information from the CTC as 
input. The resolution of the CFT is Spx/pr = 0.035(GeV/c)~^-pT- In the case of the electron 
trigger, the CFT track must be matched to a cluster in the electromagnetic calorimeter, with 
transverse energy Et > 8.0 GeV, where Et = -E- sin 6, and E is the energy of the calorimeter 
cluster. In the case of the muon trigger, the CFT track must be matched to a reconstructed 
track-segment in both sets of muon detectors. In the third level of the trigger, more stringent 
electron and muon selection criteria, which are similar to the selection criteria described in 



Section PI A| , are applied. The inclusive electron data set contains approximately 5.6 million 
events and the inclusive muon data set contains approximately 2.0 million events. These 
data are dominated by leptons from the decay of heavy fiavors {b ^ i and c ^ i) and 
hadrons that mimic the lepton signal. 
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A. Electron and Muon Identification 

Electron candidates are identified using information from botli tlie calorimeters and the 
tracking detectors. The electron calorimeter cluster in the CEM must have Et > 6 GeV. 
The longitudinal shower profile of this cluster is required to be consistent with an electron 
shower with a leakage energy from the CEM into the CHA of less than 4%. The lateral shower 
profile of the CEM cluster has to be consistent with the profile determined from test beam 
electrons. A track with pt > Q GeV/c must match the electron calorimeter cluster. This 
match is based on a comparison of the track position with the calorimeter cluster position 
determined in the CES: the difference between the extrapolated position of the track and 
the position of the cluster centroid must satisfy r|A(y9| < 1.5 cm and |Azsin^^| < 3 cm. 

To identify muons, we require a match between the extrapolated CTC track and the 
track segment in the muon chamber in both the r-ip and r-z view. The uncertainty in this 
match is taken into account and is dominated by multiple scattering in the detector material. 
The transverse muon momentum must satisfy pt > Q GeV/c. 

Finally, to ensure optimal resolution of the B hadron decay point, the electron and muon 
candidate tracks have to be reconstructed in the SVX detector. 



B. Jet Reconstruction 

Further analysis of the data sample is based on the charged particle jets in the event. 
Charged particles (instead of the more commonly used calorimeter clusters) are used to form 
jets in order to keep the electron and muon samples as similar as possible. These jets are 
found using a cone clustering algorithm. Tracks with pt > 1.0 GeV/c are used as jet seeds. 



If two seeds are within AR < 0.7 |T^, the momenta of the seeds are added together to form 
a new seed. After all possible seed merging, lower momentum tracks (0.4 < pt < 1.0 GeV/c) 
that are within AR < 0.7 of a seed are added in to form the final jets. The trigger lepton 
is always associated to a jet, and below we refer to this jet as the trigger-lepton jet. A jet 
can consist of a single track with p^ > 1 GeV/c. 



C. Secondary Vertex Reconstruction 

In order to reconstruct the time of decay in the B rest frame (the proper time), we must 
measure the point of decay with respect to the primary interaction in the lab and estimate 
the momentum of the B. Since the SVX provides only coordinates in the plane transverse 
to the beam axis, the measurement of the separation between the point of decay and the 
primary vertex is done only in the x-y plane. We refer to this separation as the decay length 
Lrcy. Only the component of the B momentum transverse to the beam axis (p^) is needed 
to calculate the proper time at decay, since the decay length is measured in the x-y plane. 

The positions of the pp interactions or "primary vertices" are distributed along the beam 
direction according to a Gaussian with a width of ~ 30 cm. In the plane transverse to the 
beam axis, these interactions follow a distribution that is a Gaussian with a width of ~ 25 fim 
in both the x and y dimensions. To reconstruct the primary event vertex, we first identify its 
2;-position using the tracks reconstructed in the VTX detector. When projected back to the 
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beam axis, these tracks determine the longitudinal location of the primary interaction with a 
precision of about 0.2 cm along the beam direction. If there is more than one reconstructed 
primary vertex in an event, the trigger lepton is associated with the primary vertex closest 
in z to the intercept of the trigger lepton with the beam line. 

The transverse position of the primary vertex is determined for each event by a weighted 
fit of all tracks with a z coordinate within 5 cm of the 2;-vertex position of the primary vertex 
associated with the trigger lepton. The tracks used in this fit are required to have been 
reconstructed in the SVX detector. First all tracks are forced to originate from a common 
vertex. The position of this vertex is constrained by the transverse beam envelope described 
above. Tracks that have large impact parameters with respect to this vertex are removed, 
and the fit is repeated. This procedure is iterated until all tracks are within the required 
impact parameter requirement. At least five tracks must be used in the determination of 
the transverse position of the primary vertex or we use the nominal beam-line position. 
The primary vertex coordinates transverse to the beam direction have an uncertainty in the 
range 10-35 /im, depending on the number of tracks and the event topology. 

The reconstruction of the B decay point (referred to below as the secondary vertex) in 
the trigger-lepton jet is based on the technique developed to identify jets formed by h quarks 
coming from t quark decay |T2|. Some modifications to this technique were necessary to 
maintain good efficiency for reconstructing the B hadron decay point in our data sample, 
since the B hadrons in this sample have substantially lower "p^ than the B hadrons from 
top quark decay. The search for a secondary vertex in the trigger-lepton jet is a two stage 
process. In both stages, tracks in the jet are selected for reconstruction of a secondary 
vertex based on the significance of their impact parameter with respect to the primary 
vertex, (io/c"do) where CTrfp is the estimate of the uncertainty on d^. The uncertainty 
includes contributions from both the primary vertex and the track parameters. The first 
stage requires at least three candidate tracks for the reconstruction of the secondary vertex. 
The trigger lepton is always included as a candidate, whether or not it satisfies the d^/adQ 
requirement. Tracks consistent with coming from the decay t[^ti~ or A" p'ir~ 

are not used as candidate tracks. Two candidate tracks are constrained to pass through 
the same space point to form a seed vertex. If at least one additional candidate track is 
consistent with intersecting this seed vertex, then the seed vertex is used as the secondary 
vertex. If the first stage is not successful in finding a secondary vertex, the second stage is 
attempted. More stringent track requirements (on do/cTdo and pt, for example) are imposed 
on the candidate tracks. All candidate tracks satisfying these stricter criteria are constrained 
to pass through the same space point to form a seed vertex. This vertex has an associated 
X^- Candidate tracks that contribute too much to the removed, and a new seed 

vertex is formed. This procedure is iterated until a seed vertex remains that has at least 
two associated tracks and an acceptable value of x^- The trigger lepton is one of the tracks 
used to determine the trigger-lepton jet secondary vertex in 96% of the events. 

The decay length of the secondary vertex Lr^y is the projection of the two-dimensional 
vector pointing from the primary vertex to the secondary vertex on the jet axis (defined 
by the sum of all the momenta of the tracks included in the jet); if the cosine of the angle 
between these two vectors is positive (negative), then L^y is positive (negative). Secondary 
vertices from the decay of B hadrons are expected to have positive L^y, while negative 
L^y vertices usually result from random combinations of mismeasured tracks. To reduce the 
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background from these false vertices, we require \Lxy/(TL^y \ > 2.0, where aL^^ is the estimated 
uncertainty on L^y. We require a secondary vertex to be associated with the trigger-lepton 
jet. This requirement leaves us with 243,800 events: 114,665 from the electron data sample 
and 129,135 from the muon data sample. The fraction of events with L^y < is 4.5% 
in the electron data sample and 5.7% in the muon-trigger sample. For reasons we discuss 
later, only events with L^y > are used in the determination of Arud and the study of the 
performance of the flavor tags. The distribution of the number of jets with total transverse 
momentum > 5 GeV/c is shown in the upper plot of Figure |l|. Approximately 60% of the 
events contain a second jet in addition to the trigger-lepton jet. The lower plot in Figure |1| 
shows the difference in azimuth between the trigger-lepton jet and the jet with the highest 
Pt in these events. A large fraction of these jets are back-to-back with the trigger-lepton jet 
as expected from the lowest-order processes that produce bb pairs. We search for secondary 
vertices in the other jets in the events as well. If an additional secondary vertex is found in 
one of these other jets, we classify this event as a "double- vert ex" event. If only the single 
secondary vertex associated with the trigger-lepton jet is found, this event is classified as 
a "single- vertex" event. The distinction between single- vertex and double- vertex events is 
important in applying the jet-charge flavor tag as described below. 



D. Determination of the B Hadron Flavor 

The next step in the analysis is to identify the flavor at production of the B hadron 
that produced the trigger lepton. We accomplish this by identifying the flavor of the other 
B hadron produced in the collision. We refer to this other B hadron as the "opposite B" in 
the text below. We first search for an additional lepton coming from the semileptonic decay 
of this opposite B. Because the pr of this lepton is not biased by the trigger, it is typically 
much smaller than the px of the trigger lepton, so we call this method of flavor identification 
the "soft-lepton tag" or SLT. The soft lepton can be either an electron or a muon. The 
lepton selection criteria are similar to the selection criteria described in Section |111 A| and 



Reference |jT2|, with additional selection criteria that use dE/dx, and pulse height in the CPR 
and CES. The soft lepton must have a track with px > 2 GeV/ c, and the invariant mass of the 
soft lepton and the trigger lepton must be greater than 5 GeV/c^. This requirement removes 
soft leptons coming from sequential semileptonic decays of charm particles produced in the 
decay of the B hadron producing the trigger lepton. The lepton identification criteria restrict 
the acceptance of the soft leptons to a pseudorapidity of \rj\ < 1.0, where r] = — ln(tan(|)). 
Approximately 5.2% of the 243,800 events contain a soft lepton candidate. 

If a soft lepton is not found, we try to identify the jet produced by the opposite B. We 
calculate a quantity called the jet charge Qjet of this jet: 

^ _ Si 9i • {Pi ■ d) .oiN 

where qi and pi are the charge and momentum of the z*^ track in the jet and d is a unit- vector 
defining the jet axis. For 6-quark jets, the sign of the jet charge is on average the same as 
the sign of the 6-quark that produced the jet, so the sign of the jet charge may be used to 
identify the flavor at production of the B hadron producing the trigger lepton. If a second jet 
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in the event other than the trigger- lepton jet has a secondary vertex, then we use this jet to 
calculate the jet charge. Double- vertex events with a jet-charge flavor tag are referred to as 
JCDV events. If only the trigger-lepton jet contains a secondary vertex, we search for a jet 
with A0 > 7r/2 with respect to the trigger lepton and pt > 5 GeV/c. If there is more than 
one jet satisfying the above criteria, we choose the jet with the highest pt- Single-vertex 
events with a jet-charge flavor tag are referred to as JCSV events. Approximately 7.5% of 
the 243,800 events above are JCDV events and approximately 42% are JCSV events. 



IV. DATA SAMPLE COMPOSITION 



The events in our selected data sample come from three sources: bb production, cc pro- 
duction, and light quark or gluon production. In each event, the trigger lepton may be a 
true lepton or it may be a hadron that mimics the experimental signature of a lepton (a 
fake lepton). The secondary vertex in the trigger-lepton jet may be a true vertex due to 
the decay of heavy flavor {b or c) or a random combination of erroneously reconstructed 
tracks that appear to form a vertex that is displaced from the primary interaction (a fake 
vertex). Light quark or gluon jets produce false vertices with L^y > with equal probability 
as false vertices with L^j^ < 0. The small fraction (<6%) of events with L^y < indicates 
that this background is small. In the analysis, we assume that all events come from heavy 
flavor [bb and cc) production. The probability of a light quark or gluon event producing 
a fake vertex and a fake lepton is negligible, although in the evaluation of the systematic 
uncertainties, we take into account the possible effects of a small amount of non-heavy flavor 
background. Below, we describe how we determine the fraction of our samples that are due 
to bb production, cc production, and fake leptons. 



A. Simulation of Heavy Flavor Production and Decay 

To understand the composition of our data, we use Monte Carlo samples of bb and cc 



production. Version 5.6 of the PYTHIA Monte Carlo generator was used to generate high- 
statistics bb and cc samples. The bb and cc pairs are generated through processes of order 
up to such as gg qq and qq qq, where q = 6 or c. Processes of order a^, such as 
gluon splitting, where gg gg is followed hj g ^ qq, are not included, but initial and final 



state gluon radiation is included. The b and c quarks are hadronized using the Peterson fl? 
fragmentation function with the parameters = 0.006 and ec = 0.06. The bottom and 
charm hadrons were decayed using version 9.1 of the CLEO Monte Carlo QQ |T^. Events 



with a lepton with px > Q GeV/c were accepted based on an efficiency parameterization |T9 
of the CFT trigger that depends on the lepton pt- The accepted events were passed through 
a simulation of the CDF detector that is based on parameterizations and simple models of 
the detector response that are functions of the particle kinematics. After the simulation of 
the CDF detector, the Monte Carlo events were treated as if they were real data. 
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B. Sources of Trigger-Electrons 



The trigger-electrons in the sample can come from three sources: heavy flavor decay 
{b —>■ e, b ^ c e, and c —>■ e), photon conversion (7 —>■ e^e~ or vr^ '~fe~^e~), or hadrons 
that fake the electron signature in the detector. The contribution from heavy flavor decay is 
discussed in section [IV D| . We attempt to identify and reject photon conversions by searching 
for the partner of the trigger-electron. We search for an oppositely charged track that forms 
a good, zero opening angle vertex with the trigger-electron. The dE/dx of this track, as 
measured in the CTC, must be consistent with the electron hypothesis. We removed 2% 
of the electron trigger sample that was identifled as photon conversions. We estimate that 
about 1% of the remaining events contain a trigger-electron from a photon conversion that 
was not identifled. To determine the fraction of events that contain a hadron that fakes 
an electron, we flt the trigger-electron dE/dx spectrum for its e, vr, K, and p content. We 
found the non-electron fraction of the sample to be (0.6 ± 0.5)%, where the uncertainty is 
statistical only. Since this background is small, we neglect it in the remainder of the analysis. 



C. Sources of Trigger-Muons 

The trigger muons in the sample can come from heavy flavor decay (6 ^ /x, 6 — > c — >■ /i, 
and c —>■ n), 71 and K decay, and from hadrons that penetrate the absorbing material in front 
of the muon chambers. The contribution from heavy flavor decay is discussed in section [iV D| . 
To study the properties of fake muon events, we used a control sample of events that only 
required a high pt CFT track in the trigger. The trigger track was treated like a trigger 
lepton, and the jet containing the trigger track was required to contain a secondary vertex. 
The L;j.y distribution of the control sample was very similar to the heavy flavor Monte Carlo 
Lxy distributions and the Lxy of the signal data samples. We conclude from this comparison 
that due to the secondary vertex requirement most of the fake muon events in the data are 
events from heavy flavor production. As described in Appendix 0, we estimate the fraction 
of events with a fake trigger muon whose dilution is zero by comparing the flavor tagging 
performance of the e and fi trigger data. The estimated fraction of events with fake muons 
that have zero dilution is (12 ± 6)%. 



D. Fraction of Data Sample due to Heavy Flavor Production and Decay 

We determine the fraction of events in the data due to bb and cc production using two 
kinematic quantities: the trigger-lepton p^^' and the invariant mass m^^ of the cluster of 
secondary vertex tracks. The quantity p^"' is deflned as the magnitude of the component 
of the trigger-lepton momentum that is perpendicular to the axis of the trigger-lepton jet. 
The trigger lepton is removed from the jet, and the jet axis is recalculated to determine p^^'- 
To calculate m'^', we assign the pion mass to all of the tracks used to form the secondary 
vertex (except the trigger lepton). We include the trigger lepton even if it is not attached 
to the secondary vertex. These kinematic quantities are effective in discriminating between 
bb and cc events because of the signiflcant mass difference between the hadrons containing 
b and c quarks (~ 3 GeV/c^). 
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Template and m distributions were obtained from the hh and cc Monte Carlo sam- 
ples. The and nf^ distributions for the data were fit to the sum of the hh and cc Monte 
Carlo templates, where the normalization for each template was a free parameter. The e 
and /i trigger data were fit separately, and the data for each trigger were divided according 
to flavor tag. The three categories were: soft-lepton (SLT), jet-charge single- vertex (JCSV), 
and jet-charge double- vertex (JCDV). The results of the Tif^ and p^^' fi^s were averaged to 
obtain the nominal values for the fraction of events from hh production The fits are 

shown in Fig. ^ and Fig. ^, and Table | gives the nominal values of F^g for the e and /i trigger 
data. The data are mostly (> 90%) from bb production. 

V. METHOD OF MEASURING THE FLAVOR TAG eD'^ AND AM^, 

As outlined in the introduction, to measure Am^ we compare the flavor of the 5° meson 
when it was produced to the flavor of the meson when it decays as a function of the 
proper decay time of the meson. 



A. Reconstruction of the B Proper Decay Time t 

To reconstruct the time of decay in the B rest frame (t), we must combine the two- 
dimensional decay length {L^y) with the component of the B momentum in the x-y plane 
(p^). The proper time is 

t = h^y:J^ (5.1) 

c ■ 

where is the mass of the B^ and c is the speed of light. The proper decay length is 
the proper time multiplied by the speed of light (ct). We do not observe all of the decay 
products of the B: the neutrino from the semileptonic decay is not detected, as well as other 
neutral decay products and charged decay products that may not have been associated with 
the secondary vertex. This means that is not known and must be estimated based on 
observed quantities and the hh Monte Carlo. 

The momentum in the transverse plane of the cluster of secondary vertex tracks p§i and 
rrf^ are the observed quantities used in the estimation of p^. The trigger lepton is included 
in the calculation of p'^ and rrf^, even if it is not attached to the secondary vertex, since we 
assume that it is a S decay product. The estimate of the B hadron transverse momentum 
p^' for an event is 

^ 1) (5.2) 

where {K) is the mean of the distribution of K = p^/p^ determined with the hh Monte 
Carlo. Figure ^ shows K distributions for two ranges of m'^^ and two ranges of p^ in the 
e-trigger hh Monte Carlo sample. For higher m'^' and p^ values, the K distribution has a 
higher mean and a narrower width: a larger fraction of the B momentum is observed so the 
observed p^ is a more precise estimate of the B momentum. To take advantage of the m'^^ 
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and p'f dependence of the K distribution, we bin the data in four ranges of m and p'f for a 
total of 16 distributions. Different sets of K distributions are used for the e and /i trigger 
data. 

Figure || shows the reconstructed proper decay length (ct) distributions for the e and 
/i trigger data. The plots of the data are compared to the expected shape from bb and 
cc production, where the bb and cc distributions were combined using in Table H and 
setting the fraction of cc events to F^c = 1 — -^^5. The lack of events near ct = is due to 
the transverse decay length significance requirement \Lxy/aL^y\ > 2. The exponential fall- 
off of the data in ct agrees with the Monte Carlo prediction. There is, however, an excess 
of events with ct < in the data. This excess is due in part to backgrounds and higher- 
order processes not included in the Monte Carlo. These backgrounds include electrons from 
photon conversions and fake muons. The higher-order processes include gluon splitting to 
bb pairs, which can produce a pair of B hadrons that are close in AR. In these events, the 
decay products of both B hadrons may be included in the same jet. This leads to secondary 
vertices that include tracks from both B hadrons, resulting in a specious measurement of 
Lxy that may be positive or negative. To verify this, we generated a sample of bb events 
using the ISAJET event generator, which includes gluon splitting. This sample showed 
an increased fraction of reconstructed vertices with negative L^y. 

B. Determination of the B Flavor at Decay 

The flavor of the B hadron associated with the trigger lepton at the time of decay is 
identified by the trigger-lepton charge, assuming the lepton is from a semileptonic B decay. 
A B hadron that contains an anti-6 quark (6) will give a positively charged lepton in a 
semileptonic decay. The trigger lepton can also originate from the decay of a charmed 
hadron produced in the decay of the B hadron, e.g. B DX,D i~X, which produces 
"wrong-sign" trigger leptons, or 5 — > J/ipX, J/ip ^ , which produces both wrong-sign 
and right-sign trigger leptons. We refer to trigger leptons from these sources as sequential 
leptons. The Monte Carlo bb sample is used to determine the fraction of trigger leptons, fseq, 
from these sequential decays. We find fseq is 9.4% in the electron sample and 13.6% in the 
muon sample. Approximately 75% of these sequential decays, i.e., 7.0% and 10.2%, produce 
wrong-sign leptons with the charge opposite to the charge from direct semileptonic decay. 
We assign a systematic uncertainty of 25% of its value to the fraction of sequential leptons, 
based on uncertainties on the branching fractions included in the CLEO Monte Carlo QQ 
and on measurements of these branching fractions at the T{4S) and the resonance. 

C. Determination of the B Flavor at Production 

To determine the flavor of the trigger-lepton B at the time of production, we attempt 
to identify the flavor of the other B in the event, and assume that the original flavor of 
the trigger-lepton B is opposite that of the other B. As described previously, we use two 
methods to obtain the flavor of the other B in the event: soft-lepton tagging (SLT) and 
jet-charge tagging (JCT). The jet-charge tag has two sub-classes: jet-charge double- vertex 
(JCDV) and jet-charge single- vertex (JCSV). The soft-lepton method is the most effective 
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(i.e., has the highest probabihty of producing a correct tag), but least efficient method. The 
jet-charge methods are less effective, but more efficient. The presence of a secondary vertex 
in the jet used for the jet charge greatly enhances its effectiveness. 



1. Quantifying the Statistical Power of the Flavor Tags 

We quantify the statistical power of the flavor tagging methods with the product eD'^, 
where e is the efficiency for applying the flavor tag, and D is the dilution, which is related 
to the probability that the tag is correct (Ptag) by 

= 2 • Ptag - 1. (5.3) 

We measure eD^ in our data. To illustrate the statistical significance of the product eD^, 
we discuss an asymmetry measurement with two types of events, a and b, where the flavor 
tagging method identifies whether the event is of type a or type b. Type a and type b 
could be "mixed" and "unmixed" decays of a neutral B meson, for example. The measured 
asymmetry ^^eas is 

where Na and Nf, are the number of events that are tagged as type a and type 6, respectively. 
The true asymmetry A is 

where A^^ and are the true number of events of type a and type b, respectively, in the 
sample. The efficiency is 

The true asymmetry is related to the measured asymmetry by 

^ = ^^measj (5-7) 

and the statistical uncertainty on the true asymmetry is 



where T is the total number of events in the sample T = Nj^ + Nj^. The statistical power of 
different flavor tagging methods varies as eD'^. 
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2. Measuring the Dilution of the Flavor Tags 



We measure the dilution of the flavor tags from our data sample. We start by deflning 
a raw dilution, -Draw^ 

^Vrs + iVws 

where N-rq and A^ws are the number of right-sign and wrong-sign events, respectively. Right- 
sign (wrong-sign) means that the charge of the trigger lepton is opposite to (the same as) 
the charge of the soft lepton or jet charge. If the charge of the trigger lepton unambiguously 
identified the fiavor at production of the B hadron, then -Draw would be equal to the true 
dilution D of the fiavor tag. All wrong-sign events would result from the fiavor tag being 
incorrect. However, since some wrong-sign events result from the trigger lepton coming 
from a B meson that mixed or from a sequential decay, -Draw is an underestimate of the 
true dilution of the fiavor tag. Backgrounds from cc production and fake leptons further 
complicate the interpretation. Nevertheless, the true dilution is approximately related to 
-Draw by a scale factor A^'^: 

D = Nd-D,^^. (5.10) 

The form of Equation p.lO| is derived in Appendix 0. We use this estimation of the true 
dilution to estimate the probability that the fiavor tag is correct on an event-by-event basis: 

Ptag = ^ (1 + AT^ ■ I^raw) , (5.11) 



This probability is used in the measurement of Am^ as described in Section |V D| . The 
dilution normalization No is determined simultaneously with Am^. 



3. Dilution for the Soft-lepton Tag 

To maximize the effectiveness of the soft-lepton flavor tag, the data are binned in the 
of the soft lepton. The p^^^ of the soft lepton is deflned in the same way as the of 
the trigger lepton. The same-sign and opposite-sign soft-lepton distributions are shown 
in Figure where the sign comparison is between the trigger-lepton charge and soft-lepton 
charge. The p^' < bin is for the case where the soft lepton is the only track in the 
jet, so that p^"^ = 0. If neither B decayed in a mixed state and if both leptons are from 
semileptonic decay, the charge of the trigger lepton would be opposite the charge of the 
soft lepton. Figure |^ shows the raw dilution -Draw for the soft-lepton tagged events. The 
raw dilution is derived from the number of same-sign and opposite-sign events in each bin 
of The dilution is lower for low because fake leptons and leptons from sequential 
semileptonic decay tend to have relatively low p^' values. The p^' dependence of -Draw was 
parameterized using the form 

i^raw(p5^^^) = A ■ (l - e-<+^) , (5.12) 
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where A and B are parameters determined from the data. The average raw dilution is used 
for isolated soft leptons that have no measurement of p^^^. The form of Equation 5.12 



IS 



empirical and was found to describe the shape of D^aw as a function of p^^^ well in the Monte 
Carlo. Parameters A and B are measured separately for the e and /i trigger data because the 
fraction of trigger leptons from sequential decay and the trigger lepton purity are different 
for the e and fi data, which affects the raw dilution. Parameters A and B are also measured 
separately for the soft-e and soft-/i tags because the fractions of fake soft leptons may be 
different for soft electrons and soft muons. The fitted values of A and B and the average 
raw dilution for isolated (no p^^^) soft leptons for soft-e and soft-/i tags in the e and fi trigger 
data are listed in Table 

The values of A and B determined from the data are assumed to describe the SLT raw 
dilution as a function of p^' for bb events. This is only an approximation since a small 
fraction (less than 10%) of the data are cc events. A cc event in which c — > i^us and 
c ^ £~us can produce opposite-sign events in which one of the leptons is associated with a 
secondary vertex, and the other lepton produces a soft-lepton tag. The soft-lepton tags in cc 
events have a much softer p^^^ spectrum than soft-lepton tags in bb events. The Monte Carlo 
predicts that cc events affect the values of A and B by an amount less than the statistical 
uncertainties on the fitted values of A and B. The effects of the approximation above are 
accounted for since the values of A and B are varied by their statistical uncertainties in 
the determination of the systematic uncertainties on Am^ and the dilution normalization 
parameters. 



4- Dilution for the Jet- charge Tag 

Figure |^ shows the jet charge distributions for single-vertex and double-vertex events. 
The data have been divided into events with a positively or negatively-charged trigger lepton 
{£). There is an anticorrelation between the sign of the jet charge and the trigger-lepton 
charge on average. The degree of separation between the £^ and i~ distributions is related 
to the raw dilution of the jet-charge flavor tag, shown in Figure ^ as a function of the 
magnitude of the jet charge |Qjet|- For double- vertex events, the presence of the second 
secondary vertex increases the probability that the jet selected for the calculation of the jet 
charge is in fact the other B in the event. This translates into a significantly higher raw 
dilution for double-vertex events. 

The IQjetl dependence of -Draw is used to predict the probability that the jet-charge tag is 
correct on an event-by-event basis, just as p^' is used for soft leptons. The IQjetl dependence 
of -Draw in the data in Figure ^ was parameterized with the form 

-Draw(IQjetl) = |Qjet| " -Dmax (5.13) 

excluding events with |(5jet| = 1 (the rightmost data point). For events with |Qjet| = 1, the 
average |Qjet| = 1 dilution is used. The slope -D^ax is determined separately for the single- 
vertex and double-vertex events in the e and fi trigger data, respectively. These slopes and 
the average raw dilution for IQjetl = 1 are listed in Table |T|. We use equation [5.11 , with 



-Draw estimated for each event using equation |5.12| for SLT events and equation [5.13| for JCT 



events, in the determination of Anid to effectively discriminate between high and low quality 
flavor tags. 
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D. Unbinned Maximum Likelihood Fit 



We use an unbinned maximum likelihood fit to simultaneously determine Am^ and the 
dilution normalization Njj. The effectiveness (eZ?^) of each flavor tag is derived from the 
measured e, -Draw? and the dilution normalization {Nij) from the fit. The fraction of 
mesons that decay in a mixed state as a function of the proper time at decay is given by 

i^mixW = ^ (1 - cos(Amd ■ t)) . (5.14) 

In a pure sample with perfect flavor tagging and proper time resolution, Fmix(i) would 
be equivalent to the fraction of same-sign events, comparing the sign of the flavors at decay 
and production. If the flavor tag is imperfect, equation p.l4| becomes 

i^mix(t) = ^tag ■ ^ (1 ~ cos(Amrf ■ t)) + Pmistag ' ^ + cos(Amd ■ t)) , (5.15) 

where Ptag (Pmistag) is the probability that the flavor tag is correct (incorrect). Using the 
relations 

Ptag = ^(l + iVz5-I)raw); (5.16) 

Pmistag ~ ~ ' ^I'aw)' (5-17) 



equation |5.15| reduces to 

Pmix(t) = i (1 - iVD ■ /^raw " COs(Amd ■ t)) . (5.18) 

Although and Arrid are correlated, the basic concept is that is determined from 
the amplitude of the oscillation in the same-sign fraction and Am^ is determined from the 
frequency. 

We determine Am^ and by minimizing the negative log-likelihood: 

- ln£ = Eln(^ss) + E HKs), (5.19) 

i j 

where Vss is the probability density for events that are tagged as same-sign, and Vqs is the 
probability density for events that are tagged as opposite-sign. Each event has three inputs 
into the likelihood: 

1. the assignment as same-sign or opposite-sign, which is based on the comparison of the 
sign of the SLT or JCT flavor tag with the charge of the trigger lepton. 



2. the estimated probability that the SLT or JCT flavor tag is correct, using equation 5.11 



with Z^raw from equation [5.12| for the SLT flavor tag or from equation |5.13| for the JCT 
flavor tag, and 



3. the decay distance L 



xy 
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In addition to the above three inputs, we use j9§i and m'^' to select the K distribution that 
is used in the determination of the reconstructed proper decay-distance. The construction 
of the probabihty densities requires several parameters. These parameters are listed in 
Table 

Both "Pss and "Pqs are the sum of several terms. First they have a term for the signal: 

T^Osi^B ) 'Ptag'^nomix ~l~ '^mistag'^mixi (5.20) 
'^Ss(-^ ) '^tag'^mix ~l~ '^mistag "^nomix i (^'^l) 



where Ptag and Pmistag are given by equations |5.16| and |5.17| , respectively, and Pnomix and 



T-'mix are given by equations |1.3| and |1.4| , respectively. Next there are terms for the other 
B hadrons, including terms for both B^ and h baryons, as well as a term for 5°. The terms 
for B^ are 

^os(S+) = Ptag— e-*/^^+; (5.22) 

Pss(5+) = Pmi.tag— e-*/^^+ ; (5.23) 
'Tb+ 

where tb+ is the lifetime of the B^ , and t is the proper decay-time. The terms for h baryons 
are similar, except that tb+ is replaced by the lifetime of the h baryons, Tbaryon- The terms for 
i?^ are similar to the terms for B^ , except that Am^ is replaced by Am^, which is assumed 
to be very large {i.e. beyond our experimental sensitivity) so that these terms effectively 
look like 

Vo,{Bl) = Pss(i?°) = ^e-'^^-2. (5.24) 

2tbo 

The values of the lifetimes of the B hadrons and the value of Arris used in the probability 
densities are listed in Table |IV|. Each term for the various B hadrons is multiplied by the 
expected relative contribution, /^o, /bo, and /baryon, of these various hadrons to the 
data sample. The production fractions are renormalized to take into account the different 
semileptonic branching fractions of the various B hadrons: the semileptonic widths are 
assumed to be identical, so the semileptonic branching fractions are scaled to agree with the 
relative lifetimes. 

In addition to the terms for direct semileptonic decay, there are terms that take into 
account the contribution from sequential decays. The fraction of sequential decays fseq is 
listed in Table 75% of these sequential decays produce leptons with a sign opposite to 
direct semileptonic decay. 

The terms that take into account the contribution from cc events are similar to the terms 
for B^, except that we use an assumed flavor tagging dilution and an effective lifetime Tcc- 
The true flavor tagging dilution is not known a priori for both bottom and charm decays. 
Just as for B decays, we do not consider the Monte Carlo reliable for predicting the JCT 
or SLT dilution for charm decays, so we use assumed values for the charm dilution and 
varied them by the maximum possible amount for the systematic uncertainties. For the 
JCT, we expect the dilution for charm decays to be worse than for bottom decays, therefore 
we assume Dcc/Df^i = 0.5 for the JCT, and vary Dcc/Df^i from to 1 in the evaluation of 
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the systematic errors. The ratio Dec/ Df,i is used to rescale the predicted bb dilution for the 
event, based on \Qjet\ using Equation |5.13| , to give the predicted cc dilution for the event. 
The SLT dilution for charm decays could be anything from to 1 depending on the fraction 
of fake soft leptons in cc events. Unlike bottom decays, the SLT dilution for charm decays 
does not fall near p^^' = due to soft leptons from sequential decays, therefore we assume 
Dec = 0.5, independent of for the SLT, and vary Dcc/Di^i from to 1 in the evaluation 
of the systematic errors. 

The effective lifetime Tcc = 1.53 ps is determined from the cc Monte Carlo samples, 
where the proper time at decay was reconstructed using the K distributions from the bb 
Monte Carlo and equation [5.1| . The relative fractions of hb {Fjji) and cc (-Fee) ? which can be 
determined from Table |, multiply the terms in the likelihood corresponding to bb and cc 
production, respectively. Finally, for the case of the //-trigger data, a term for fake muons 
was included with the lifetime and zero dilution. The relative fraction of fake muons 
-^fake fi (discussed in Appendix ^ is listed in Table |V[ In this case, and Fee are scaled 

by 1 - i^f°ke 

The probability densities Pss(^) and Pos(^) are functions of the true proper time at decay 
t. We take into account the experimental resolution on t by convoluting Pss{t) and Pos(^) 
with Lrcy and p^' resolution functions. Distributions of SL^y = L^-j^ (measured) — Lxy{tTue) 
from the e and fi trigger bb Monte Carlo samples are parameterized with the sum of three 
Gaussians, a SL^y > exponential, and a SL^y < exponential. The SL^y distributions 
and their parameterizations are shown in Figure |T0| for the e and fi trigger data. The 
K distributions, like those shown in Figure are used as p^' resolution functions. The 
Lxy and p^' resolution functions describe the t smearing for bb events from direct production 
{gg — i> bb and qq bb) . They do not describe backgrounds such as photon conversion events 
(for the e trigger) or events with a fake trigger muon. They also do not describe bb events 
from gluon splitting {gg —>■ gg followed hj g bb), which tend to have worse t resolution 
due to secondary vertices that include decay products from both B hadrons. While the 
total amounts of these backgrounds are reasonably small, they become more important near 
ct = and for very large ct values (beyond several B lifetimes). For these reasons, we only 
use events with 0.02 cm < ct < 0.30 cm in the fit for Am^ and the dilution normalization 
Nd- 



VI. FIT RESULTS 

The free parameters in the unbinned maximum likelihood fit are the mass difference Am^ 
and the dilution normalization factors No- There are six dilution normalization factors: 

1. jQ^y and -/V/^ jcsv- t^e dilution normalization factors for the jet-charge flavor tag 
in the case that the jet does not contain a secondary vertex for the electron-trigger 
and muon-trigger data samples, respectively. 

2. and N'^ j^j^y-. the dilution normalization factors for the jet-charge flavor tag 
in the case that the jet does contain a secondary vertex for the electron-trigger and 
muon-trigger data samples, respectively. 
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3. Nf) gLrp and A^^^^ gL^: the dilution normalization factors for the soft-lepton flavor tag for 
the electron-trigger and muon-trigger data samples, respectively. 

To determine the dilution normalization factors needed to calculate the flavor tag eD^ values 
(see Section |V1 C|) , the dat grouped into four subsamples: 



1. e-trigger, soft-lepton flavor tag, 

2. e-trigger, jet-charge flavor tag, 

3. /i-trigger, soft-lepton flavor tag, and 

4. /i-trigger, jet-charge flavor tag. 

There is some overlap of events in these four subsamples as some events have both a SLT 
and a JCT. About 20% of the events with soft-lepton tags did not pass the single-lepton 
Level 2 trigger and came, instead, from a dilepton trigger. In order for the SLT efficiency 
to be well defined, we require that the SLT events pass the single-lepton trigger in data 
subsamples 1 and 3, which are used to determine Nf^^^^r^, ^d,sijt^ ^^"^ the SLT efficiencies 
needed for calculating the SLT eD"^. 

The fit results for the individual flavor taggers are given in Table |V[ The dilution 
normalization factors are given for two cases: Am^ free to float in the fit and Am^ fixed to 
the world average |]2T]] ( 0.474 Ti ps~^). The value of Am^ is held fixed to the world average 
because Am,^ and the No factors are correlated. The correlation coefficients between Am^ 
and the No constants range from 0.55 to 0.81. Fixing Am^ reduces the statistical uncertainty 
on the dilution normalization constants and removes any bias from statistical fluctuations 
that pull Am^ high or low. The No factors determined with Am^i fixed to the world average 
are used in the calculation of the flavor tag eD^ values. The A^^^ factors in Table |V| are 
consistent with our expectations from the composition of the data (see Appendix 

To determine Arrid, we fit all four data subsamples simultaneously. For events with both 
a SLT and a JCT, we use the SLT because it has significantly higher average dilution. The 
results of the simultaneous fit of the e and /i trigger data using both flavor tagging methods 
are listed in Table [V^. We find Am^^ = 0.500 ± 0.052 h ps~^, where the uncertainty is 
statistical only. This is consistent with the world average value (0.464 ± 0.018 h ps^^) 0. 
Using the SLT tag in doubly flavor tagged events removes events with higher-than-average 
dilution from the JCT. This results in lower values of No for the jet-charge flavor tag than 
found in the fits of the individual samples listed in Table |V|. 

Figure 11 shows the fraction of same-sign events as a function of the reconstructed proper 
decay length ct. The points with error bars are the data. The curve is a representation of 
the results of the fit for Am^ and No- Figure has been included to illustrate the clear 
evidence of mixing in the data. However, it does not contain all of the information that 



goes into the unbinned likelihood fit. In Figure 11, all events are treated equally. In the fit. 



events are effectively weighted based on their estimated dilution. 



A. Check of Fitting Procedure 

To check our fitting procedure, we used a fast Monte Carlo which generated hundreds 
of data samples, each with the same statistics, tagging dilution, and t resolution as the 
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real data. Figure shows the fit results of 400 fast Monte Carlo samples, representing 
the SLT flavor tagged, e-trigger data. The top row of plots show distributions of the fitted 
values of Am^ and Nj:, for the 400 samples. The arrows indicate the values of Am^ and 
A^D with which the samples were generated. The mean of each distribution is consistent 
with the generation value. The middle row of plots show the distributions of the statistical 
uncertainty on Am^ and No- The arrows indicate the statistical uncertainty from the fit to 
the SLT flavor tagged, e-trigger data. The statistical uncertainties on Am^ and No for the 
data are near the most probable values from the fast Monte Carlo samples. The bottom row 
of plots show the distributions of the deviation of the fitted value from the generation value 
divided by the statistical uncertainty for Am^^ and N^. These distributions have a mean of 
zero and unit width, which confirms that Am^ and Na are unbiased and that the statistical 
uncertainty is correct. 



B. Systematic Uncertainties 

To determine the systematic uncertainty on Am^, the fixed input parameters in Table ^ 
were varied by the amounts listed in this Table. To determine the systematic uncertainty 
associated with the dilution parameterizations, the parameters describing the -Draw depen- 
dence on IQjetl and of the soft lepton were varied by their statistical uncertainties. A 
systematic uncertainty to account for the possible presence of fake secondary vertices from 
non-heavy flavor backgrounds as well as heavy flavor events (gluon splitting) was determined 
using a combination of data and Monte Carlo. The observed excess in the data over the 
combined contribution of hb and cc (see Figure ^ was used to define the shape of this back- 
ground. This shape was included in the likelihood function with a dilution that varied from 
to D,i. 

The systematic uncertainty assigned to the variation of each input parameter was the 
shift in Amrf from the fit result with the nominal values of the input parameters. The 
total systematic uncertainty is the sum in quadrature of these shifts. The same procedure 
was applied to the four data subsamples (e SLT, e JCT, /i SLT, /i JCT) to determine the 
systematic uncertainties for the dilution normalization factors. 

The procedure described above was checked for the largest individual contributions to the 
systematic uncertainty using fast Monte Carlo samples. Samples generated with a variation 
on one of the input parameters (e.g. tb+/tbo = 1.02 + 0.05) were fit using the nominal 
parameter (e.g. tb+/tbo = 1.02). The average bias on the fitted values for the fast Monte 
Carlo samples was consistent with the deviation on the fitted values observed when the data 
are fit with a fixed parameter variation. 



Table [VI I| lists the individual contributions to the systematic uncertainty on Am^ 
(±0.043 h ps^^). The largest single contribution (±0.032 h ps^^) is the unknown soft- 
lepton flavor tag dilution for cc events. The SLT cc dilution was varied over the full possible 
range (0 to 1). If the assumed cc dilution in the fit is lower (higher) than its true value, 
the fraction of same-sign events at small t will appear lower (higher), which biases Arrid low 
(high). The second largest contribution (±0.021 h ps~^) is the uncertainty on the lifetime 
ratio of the and the mesons. If the value of tb+ / tbo in the fit is higher (lower) than 
its true value, the fraction of same-sign events at short t values will appear larger (smaller), 
which biases Am^ high (low). The third largest contribution (±0.009 h ps~^) is the uncer- 



23 



tainty on the SLT raw dilution parameterization: Draw as a function of soft-lepton p^^'. The 
size of the variation for the SLT is shown by the dashed curves in Figure |^. The raw dilution 
for events with no p^"' nieasurement was independently varied by the statistical uncertainty 
on the raw dilution for no-p^?' events. 

Tables |V111| and |I^ list the contributions to the systematic uncertainty on the dilution 
normalization factors for the e and ^ trigger data respectively. As for the systematic un- 
certainty on Am^, the largest contribution comes from the unknown dilution for cc events. 
The fraction of Bg events, which we assume are half same-sign and half-opposite sign, and 
the fraction of events where the trigger lepton is from a sequential decay both affect the as- 
signment of same-sign events. If these fractions are low (high), more (less) same-sign events 
will be attributed to mistags, thus they are strongly coupled to the dilution normalization. 
The uncertainty on the raw dilution parameterizations also has a large effect on the dilution 
normalization systematic uncertainty. 



C. Flavor Tag eD"^ 

The measurement of the statistical power tD^ of the flavor tagging methods is done 
in two steps. First, the raw tD^ is calculated using the raw dilution rather than the true 
dilution. Then, the raw tD^ is rescaled by N]^, which translates the raw dilution to the true 
dilution. The Nr, factor for each flavor tagging method in the e and /i samples is determined 
in the unbinned maximum likelihood fit with Am^ fixed to the world average ||2T| (0.474 
h ps~^). These Nd values are given in Table 0, where the first uncertainty is statistical and 
the second systematic. 

The raw or uncorrected value of eD^ is obtained by summing the e-Draw values in the 
bins of either p^^^ for the SLT or \Q]et\ for the JCT. The efficiency for each bin is defined as 
the number of events in the bin divided by the total number of events before flavor tagging. 
Table ^ lists the total efficiency, raw tD^, dilution normalization, and true eD^ for each 
of the flavor tagging methods. Taking the average of the e and /i trigger data, we find 
eD"^ to be (0.78 ± 0.12 ± 0.08)% for the JCT and (0.91 ± 0.10 ± 0.11)% for the SLT where 
the first uncertainty is statistical and the second systematic. These eD^ values are about 
one order of magnitude lower than typical flavor tagging techniques employed on the Z° 
resonance [0], however the large hb cross section in pp collisions at ^/s = 1.8 TeV yields hh 
samples that are about one order of magnitude larger at CDF than those collected on the 

resonance. The statistical uncertainty for our measurement of Am^ is still competitive 
with similar measurements on the resonance, since our smaller eD^ is compensated by 
our larger sample size. 

The values of eD^ for these flavor tags depend on the data sample in which they are used. 
In particular, during next run of the Tevatron, we will collect large samples of / — > 
J/ipRg (for the precise measurement of the CP asymmetry parameter sin(2/3)) and hadronic 
-Bg decays (for the precise determination of Arris). The triggers used to collect these data 
samples will be different from the inclusive lepton trigger used to collect the data for this 
analysis. As a result, the B hadron production properties [e.g., px of the B) are different, 
and this affects eD^. Despite these differences, the results in this paper demonstrate that 
both the jet-charge and the soft-lepton flavor tagging methods are viable in the environment 
of pp collisions. 



24 



VII. SUMMARY 



We have measured Am^ using soft-lepton and jet-charge flavor tagging methods. This 
is the first apphcation of jet-charge flavor tagging in a hadron-collider environment. The 
flavor at decay was inferred from the charge of the trigger lepton, which was assumed to be 
the product of semileptonic B decay. The initial flavor was inferred from the other B in the 
event, either using the charge of a soft lepton or the jet charge of the other B. The proper- 
time at decay for each event was determined from a partial reconstruction of the decay 
vertex of the B that produced the trigger lepton and an estimate of the B momentum. The 
value of Arriii was determined with an unbinned maximum likelihood flt of the same-sign 
and opposite-sign proper-time distributions (comparing the sign of the trigger-lepton charge 
and the flavor tag). The statistical power of the flavor tagging methods was measured in 
the unbinned maximum likelihood flt by fltting for a scale factor Nd, for each of the flavor 
tagging methods, which is the ratio of the raw dilution and the true dilution. 

We find Airid = 0.500 ± 0.052 ± 0.043 h ps^^, where the first uncertainty is statistical 
and the second systematic. This is consistent with the world average value of 0.464 ± 
0.018 h ps~^ [0 and competitive in precision with other individual measurements of Am^. 
We quantify the statistical power of the fiavor tagging methods with eD^, which is the tagging 
efficiency multiplied by the square of the dilution. We find eD'^ to be (0.78 ± 0.12 ± 0.08)% 
for the jet-charge fiavor tag and (0.91 ± 0.10 ± 0.11)% for the soft-lepton fiavor tag, where 
the first uncertainty is statistical and the second systematic. These eD^ are much lower than 
what has been achieved in experiments on the Z° resonance, however we have demonstrated 
that the much higher bb cross-section at the Tevatron {pp, ^/s = 1.8 TeV) can be used 
to compensate for the disadvantage in eD^. The jet-charge and soft-lepton fiavor tagging 
techniques will be important tools in the study of CP violation in the up-coming run of the 
Tevatron. 
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APPENDIX A: THE DILUTION NORMALIZATION FACTOR 



The true dilution D of a flavor tagging method is defined as 



D = 2-R 



tag 



(Al) 



where Ptag is the probabihty that the flavor tag is correct. An equivalent expression for D 
is 



D 



M 



Nt + N, 



M 



(A2) 



where Nt (Nm) is the number of correct (incorrect) tags in a sample of A^totai 
events. The raw dilution is deflned as 



M 



D 



OS 



ss 



- Nos + Ns ' 



(A3) 



where Nqs (Nss) is the number of opposite-sign (same-sign) events in the sample, comparing 
the trigger-lepton charge with either the sign of the soft-lepton charge or the jet charge. If 
the data were pure bb with no B mixing and all of the trigger leptons were from direct B 
decay, all opposite (same) sign events would be correctly (incorrectly) flavor-tagged. That 
is, we would have Nt = Nqs, Nm = Nss, and D = -Draw There are, however, several things 
in the data that break the Njn = Nqs and Nm = Nss assumptions. They are 

• B mixing: If the trigger lepton is from a B hadron that decays in a state opposite 
its original flavor, the trigger-lepton charge will have the "wrong" sign. In this case, 
events with the correct flavor tag are same- sign. 

• Sequential decays: The charge of trigger leptons from sequential B decay 
{b ^ c ^ £ s X) is opposite that of direct B decay. For trigger leptons that are from 
sequential decay, events with the correct flavor tag are same-sign, if the trigger-lepton 
B did not mix. 

• cc events: Events from cc production may have a non-zero dilution that is not the 
same as the dilution from bb events. 

• Fake leptons: The e-trigger data has essentially no fake trigger electrons. However, 
about 12% of the /i-trigger data have a hadron that faked a muon, whose charge is 
random (see Appendix |Bp . 

If there were no fake leptons, the number of opposite-sign and same-sign events from bb 
production are given by 
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where is the fraction of trigger leptons in bb events that are from sequential decay in 
which the trigger-lepton charge has the "wrong" sign, x' is the effective^] probabihty that 
the B hadron that produced the trigger lepton decayed in a mixed state, and Nq^ (N^'g) is 
the number of same-sign (opposite-sign) bb events. For events from cc production we have 



T^CC 
T^CC 



T^CC 



(A6) 
(A7) 



where Nq^ {Nf-g) is the number of same-sign (opposite-sign) cc events. Using Equations |A3| , 
ISE, RE, and Wj\ the raw dilution can be written as 



(1 - 2x0 (1 - vzi)m' - m^) + m - n\ 



(A8) 



M 



If a fraction of the events (-Fffke have a fake trigger lepton whose charge-sign is random, 
these events will have a raw dilution of zero since the number of same-sign fake-lepton events 
will equal the number of opposite-sign fake-lepton events. Taking fake leptons into account 



gives 



-^fakc 



(A9) 



M 



Using Equation [X2| , we define the true flavor tagging dilution in bb and cc events as 

- (AlO) 



and 



N^' + AT^ 



NS? + A^^f 



(All) 



respectively. The fraction of events from bb and cc production are defined by 



Fbb = 



M 



and 



+ A^; 



" N^^ + Nf, + N^^ + AT-" 



(A12) 



(A13) 



respectively. Combining Equations [A8| , [AlO] , |A11| , [A12| , [A13| gives 



^ It is an "effective" probability because our secondary vertexing method is inefficient for low 
values of t, which causes x' to be larger than x 0- 
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where we have defined the dilution normahzation factor as 



— (1 -^fakc l) 
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(A15) 



Equation |A15| can be used to calculate the expected values for the A^^^ parameters. For this 
calculation, we will assume 



x' ~ 0.20 from the Monte Carlo. 

/s'^q(e-trigger) = 0.07 and /,^^(/x-trigger) = 0.10 using f^l = 0.75 x /,eq and the 



values in Table IV. 



The Ff^i values are given in Table |. We also use F^ 



66- 



• For the JCT, we assume Dcc/Df^i^JCT) = 0.5. Using the average SLT dilution and 
the assumption that Dcc(SLT) = 0.5, we estimate Dec/ Df,i{SLT) ~ 1-3. 

Using the numbers above, we find 

• A^B.SLT ~ 1-8. 
. Ar^,sLT^2.1. 

a Kfe Tvre ~ 1 Q 



-^D,JCSV ~ ^D.JCDV ~ 2.4. 
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APPENDIX B: FRACTION OF FAKE TRIGGER MUONS 



As is stated in Section [IV (J| , we believe that most of the fake trigger-muons in the data 



are from heavy flavor decay. There may be some correlation on average between the sign 
fake muon charge and the B flavor at decay, however we assume that this correlation is 
smaller than that of real trigger-muons. Fake muon events are divided into two groups: 

1. Fake muon events whose dilution is the same as the dilution for real muons. 

2. Fake muon events whose dilution is zero. 

We treat group 1 as if they are real trigger-muon events. We treat group 2 as if they are bb 
events with a flavor tagging dilution of zero. 

We determine the fraction of events with fake muons that have zero dilution in the /x- 
trigger data by assuming that the true flavor tagging dilution is the same in the e and n 
trigger data. The raw dilution (Draw), which assumes all opposite-sign (same-sign) events 
are tags (mistags), is different for the e and /i trigger data for the following reasons: 

1. The fraction of real trigger- leptons in bb events that are not from direct b ^ i decay 
(/seq) is 9.4% for the e-trigger data and 13.6% for the /i-trigger data. 

2. The fraction of events from cc production F^c is slightly different (see Table |ip. 

3. We estimate that only 1% of the trigger electrons are fake, while, as shown below, 
about 10% of the /i-trigger events contain a fake muon. 

We can correct Z^raw for (1) and (2) using the equation 

/T-)/ \ _ (-Draw) 

^ ''"^^ F,i-{l-2fZ\ + F,,.D,,/D,j,) ^ ^ 

where f^^^ is the fraction of non b ^ i decays that have the "wrong" sign. The Monte 
Carlo gives fZl = 0.75 ■ /seq. The values of {D'^^J for the SLT, JCSV, and JCDV flavor 
tagging methods in the e and /i-trigger data are given in Table ^11 The weighted average 



of {D[.^M) KD',^M) for the SLT, JCSV, and JCDV flavor tagging methods gives 1.14 ± 
0.08. The fraction of events with fake muons that have zero dilution can be extracted using 

Equation ^ gives -Ffake^ = 12 ± 6 %. The relatively large uncertainty on -Ffakc^ gives a 
signiflcant systematic uncertainty on the dilution normalization iV^) for the flavor tags (see 



Table p(| ), however, the contribution to the systematic uncertainty on Am^ is relatively 
small. 
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TABLES 

TABLE I. The fraction of events from bb production Fj,^ for the e and /j, trigger data. The first 

two cohimns give the result for from fitting the p^^ and m^' spectra to a Uncar combination of 
bb and cc templates. The error in the first two columns is statistical only. The last column gives 
the average of the p^^^ and m*^' results. The error on the average takes into account the difference 
in the p^^ and m^' fit results as a systematic error. The fraction of events from cc production is 
given by Fee = 1 - F^l- 



Sample 


p'rf Fit 


rrf^ Fit 


Average Value 


JCSV(e) 


89.6 ± 0.6 % 


92.2 ± 0.5 % 


90.9 ± 1.3 % 


JCDV(e) 


94.5 ± 1.5 % 


96.7 ± 1.0 % 


95.6 ± 1.5 % 


SLT(e) 


91.3 ± 1.7 % 


94.4 ± 1.5 % 


92.9 ± 1.7 % 


JCSV(ai) 


90.1 ± 0.5 % 


89.0 ± 0.5 % 


89.6 ± 1.3 % 


JCDV(Ai) 


98.7 ± 1.3 % 


97.6 ± 1.2 % 


98.2 ± 1.3 % 


SLT(/x) 


92.8 ± 1.6 % 


91.5 ± 1.4 % 


92.2 ± 1.6 % 



TABLE IL Fitted values for the SLT raw dilution parameterization constants {A and B) and 
the raw dilution for isolated soft leptons (no-p^') for soft e and pL tags in the e and jJL trigger data. 
The errors are statistical. The SLT raw dilution is parameterized as a function of the soft lepton 
p^rf^ with the functional form -Draw(PT') = A • ^1 — 6"^^'+^^ where p'rf^ is in GeV/c. 





e-Trig 


ger 


/i- Trigger 






Soli-e 


Sol'l-// 


Sofl-f 




A 


0.41 ± 0.05 


0.46 ± 0.03 


0.42 ± 0.05 


0.36 ± 0.03 


B 


0.31 ± 0.13 


0.43 ± 0.06 


0.35 ± 0.11 


0.25 ± 0.09 


no-p^^ 


0.36 ± 0.04 


0.25 ± 0.03 


0.18 ± 0.04 


0.18 ± 0.03 



TABLE IIL Constants for the jet-charge tag raw dilution parameterization as a function of 
IQjetl- The parameterization has the form -Draw = -Dmax • IQjetl except for events where IQjetl = !> 
in which case the average -Draw is used. 





e-Trigf 


^er 


/x-Trigj 


^er 




JCSV 


JCDV 


JCSV 


JCDV 


-Dmax 


0.083 ± 0.010 


0.34 ± 0.02 


0.060 ± 0.009 


0.29 ± 0.02 


-Draw) IQjetl ~ 1 


0.091 ± 0.014 


0.18 ± 0.03 


0.074 ± 0.014 


0.12 ± 0.03 
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TABLE IV. The fixed input parameters used in the fit. The first column hsts the parameters, 
which are described below. The second column lists the value of the parameter assumed for the 
fit, with an error on this parameter, if appropriate. The third column lists the variation of the 
parameter for the purposes of evaluating the systematic errors on Am^ and Nd; itlo" means that 
the parameter was varied by the error listed in the second column. Finally the fourth column 
lists the source of the parameter. The fixed input parameters are the B hadron lifetimes (rj), the 

mass difference (Attt-s), the relative B hadron production fractions (/j), the fraction of trigger 
leptons from bb sources other than direct B decay (fseq), the flavor tag dilution value for cc events 
(Dec) , an effective lifetime for cc events (tcc) , and the fraction of //-trigger events where the trigger 

fake n ) ' 



muon is really a hadron and has a dilution of zero {F^ 



Parameter 



Value 



Variation 



Source 



tb+/tb" 



1.56 ±0.06 ps 

1.02 ±0.05 

1 61+°-^° ns 
-'■■0-'--o.09 

1.14 ±0.08 ps 



±lcr 
±lo- 

±1CT 



|21] 



m 

m 
12111 



Am, 



700 ?ips 



-1 



down to 6.7 ?ips 



-1 



Assumption 



Ibo 



37.8 ± 2.2 % 
37.8 ± 2.2 % 

13.2 ±4.1 % 



±lo- 

±1(7 
±1(7 
±1(7 



±0.25 X /seq 
±0.25 X /seq 



m 



m 
n 



seq 



e: 9.4 % 
//: 13.6 % 



Monte Carlo 
Monte Carlo 



JCT D,,/D,-, 
SLT D,s 

'Tcc 



0.5 
0.5 

1.53 ±0.20 ps 



to 1 
to 1 

±1(7 



Assumption 
Assumption 
fit to cc MC 



7?o 

fake fi 



12 ± 6 % 



±1(7 



e/fj, Dj-i^^ comparison 



TABLE V. Results of the maximum likelihood fit for Anid and the flavor tag dilution nor- 
malization factors for the individual flavor taggers separately in the e and fi trigger data. The 
first error is statistical and the second is systematic. The evaluation of the systematic errors is 
discussed later in the text. The dilution normalization factors are given for two cases: Anid free 
to float in the fit and Anid fixed to the '96 world average [18] ( 0.474 h ps^^). 

SLT JCSV JCDV 



e-trigger 

Nd, Ama free 1.66 ± 0.12 ± 0.18 

Nd, Amd fixed 1.72 ± 0.08 ± 0.11 

Anid {h ps~^) 0.45 ± 0.08 ± 0.05 

/i-trigger 

Nd, Amd free 2.05 ± 0.19 ± 0.37 

Nd, Amd fixed 2.01 ± 0.13 ± 0.22 

Amd (?i ps"^) 0.50 ± 0.09 ± 0.05 



1.84 ± 0.23 ± 0.19 1.71 ± 0.16 ± 0.11 

1.88 ± 0.20 ± 0.15 1.76 ± 0.13 ± 0.09 

0.42 ± 0.09 ± 0.03 



2.86 ± 0.40 ± 0.43 2.52 ± 0.28 ± 0.25 

2.41 ± 0.29 ± 0.39 2.14 ± 0.33 ± 0.25 

0.68 ± 0.11 ± 0.04 
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TABLE VI. Results of the maximum likelihood fit for Am^ and the flavor tag dilution normal- 
ization parameters {Njj) for the e and fj, trigger data using both soft-lepton and jet-charge flavor 
tagging. The first error is the statistical error on the parameter, determined by the fit, and the 
second error is the systematic error, which is determined using the prescription described in the 
text. The dilution normalization parameters (Nd) for the JCT are lower than those in Table ^ 
because the SLT is used for double flavor tagged events (SLT and JCT), which lowers the average 
dilution for the JCT. Only the statistical errors are given for the dilution normalization parameters. 

Parameter Fit Result 

Anid jh ps"^) 0.500 ± 0.052 ± 0.043 

iV|, jcsv 1.43 ±0.22 

^^AJCDV 1.53 ±0.16 

A^I^SLT 1.72 ±0.10 

K,jcsv 1.99 ±0.32 

<,JCDV 2.00 ±0.22 

% SLT 2.05 ±0.16 
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TABLE VII. The individual contributions to the systematic error on Am^. The value of SAnid 
is the amount by which the fitted value of Am^ shifted when a parameter was varied. The size 



of the variations are given in Table |V| and in the text. The parameters that were varied are: 
the B hadron lifetimes (rj), the assumed value of Atti-s, the and A;, production fractions (fi), 
the fraction of trigger- leptons from sequential decay (/seq), the assumed dilution for cc events, the 
effective lifetime for cc events (tcc) , the fraction of trigger muons that are fake and have zero dilution 
(-^fake /i)' normalization for the fake- vertex component, the assumed dilution for fake- vertex 
events [D-py / D^i), the fraction of cc events Fee, and the parameterization of Draw The total 
systematic error is the sum of the individual contributions in quadrature. 



Parameter 


dArrid [ri ps ) 




0.004 


1 

'TB+n BO 


0.021 


Tbo 


0.005 




0.005 




r\ r\r\ a 

0.004 




0.001 




0.007 


/seq 


0.004 


JCT De-clD,i 


0.002 


SLT Dec 


0.032 


Tcc 


0.006 


fake fj, 


0.001 


Fake- vertex component 


0.006 




0.003 


Fee, e-trigger 


0.006 


Fee, /U-trigger 


0.003 


JCT Draw parameterization 


0.003 


SLT -Draw parameterization 


0.009 


Total systematic error 


0.043 
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TABLE VIII. The individual contributions to the systematic error on the dilution normaliza- 
tion parameters for the e-trigger data Nf^ •. The value of SNf^ • is the amount by which the fitted 



value of Nf) • shifted when a parameter was varied. The size of the variations are given in Table |V 
and in the text. These systematic errors are for the case where Am^ is allowed to float in the fit. 
The parameters that were varied are: the B hadron lifetimes (rj), the assumed value of Amg, the 
and production fractions (/,), the fraction of trigger-leptons from sequential decay (/seq), 
the assumed dilution for cc events, the effective lifetime for cc events (tcc), the normalization for 
the fake- vertex component, the assumed dilution for fake- vertex events {Dpy /Dj^i), the fraction 
of cc events Fee, and the parameterization of Draw The total systematic error is the sum of the 
individual contributions in quadrature. 



Parameter 








T go 


0.01 


0.01 


0.01 


TB+frB^ 


0.01 


0.01 


0.01 




0.01 


0.01 


0.01 




0.00 


0.00 


0.00 


Anis 


0.02 


0.01 


0.01 


fso 


0.05 


0.05 


0.05 




0.01 


0.01 


0.01 


/seq 


0.06 


0.06 


0.06 


Assumed cc dilution 


0.12 


0.05 


0.15 


Tec 


0.00 


0.00 


0.01 


Fake- vertex component 


0.02 


0.01 


0.02 




0.06 


0.03 


0.04 


F - 
^ cc 


0.00 


0.01 


0.01 


Dram parameterization 


0.10 


0.05 


0.06 


Total Systematic Error 


0.19 


0.11 


0.18 
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TABLE IX. The individual contributions to the systematic error on the dilution normalization 
parameters for the /x-trigger data Nj^ ■ . The value of SN'^ ■ is the amount by which the fitted value 



of • shifted when a parameter was varied. The size of the variations are given in Table IV 



and in the text. These systematic errors are for the case where Am^^ is allowed to float in the 
fit. The parameters that were varied are: the B hadron lifetimes (rj), the assumed value of Am^, 
the and A;, production fractions (/j), the fraction of trigger-leptons from sequential decay 
(/seq)) the assumed dilution for cc events, the effective lifetime for cc events (tcc), the fraction 
of trigger muons that are fake and have a dilution of zero (-^f^ke the normalization for the 
fake-vertex component, the assumed dilution for fake- vertex events {Dpy / Dy^), the fraction of 
cc events Fee, and the parameterization of -Draw The total systematic error is the sum of the 
individual contributions in quadrature. 



r^arameter 




"^^D,JCDV 




r go 


0.00 


0.00 


0.01 


TB+frB^ 


0.02 


0.02 


0.00 




0.01 


0.01 


0.01 




0.02 


0.02 


0.01 


Arris 


0.02 


0.01 


0.01 


/so 


0.08 


0.07 


0.06 




0.02 


0.01 


0.01 


/seq 


0.15 


0.14 


0.13 


Assumed cc dilution 


0.28 


0.03 


0.26 


'Tec 


0.01 


0.00 


0.01 


i?o 

fake /i 


0.19 


0.17 


0.15 


Fake vertex component 


0.00 


0.01 


0.01 




0.13 


0.08 


0.06 


F - 

^ ee 


0.00 


0.01 


0.02 


Draw parameterization 


0.17 


0.03 


0.07 


Total Systematic Error 


0.43 


0.25 


0.37 
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TABLE X. The statistical power eD^ for the flavor tagging methods used: Jet-Charge Single 
Vertex (JCSV), Jet-Charge Double Vertex (JCDV), and Soft-Lepton Tag (SLT). Results for the 
e and /i trigger data are shown in separate rows. The sum is over bins of for the soft-lepton 
data and |Qjet| for the jet-charge data, as shown in Figures and respectively. The square of 
the dilution normalization factor Nd is used to rescale the J2i ^i-^raw i value to give J2i ^i^f- The 
first error is statistical, the second systematic. 



Sample 


Total e 


T e D"^ 

i—ii ^* raw I 






JCSV (e) 
JCDV (e) 
SLT (e) 


41.55 ± 0.14 % 
7.44 ± 0.08 % 
4.38 ± 0.06 % 


0.077 ± 0.016 % 
0.159 ± 0.023 % 
0.329 ± 0.033 % 


1.88 ± 0.20 ± 0.15 
1.76 ± 0.13 ± 0.09 
1.72 ± 0.08 ± 0.11 


0.27 ± 0.06 ± 0.04 % 
0.49 ± 0.10 ± 0.05 % 
0.97 ± 0.13 ± 0.12 % 


JCSV ifi) 
JCDV in) 
SLT (/i) 


43.81 ± 0.14 % 
7.66 ± 0.07 % 
4.54 ± 0.06 % 


0.048 ± 0.012 % 
0.113 ± 0.018 % 
0.210 ± 0.026 % 


2.41 ± 0.29 ± 0.39 
2.14 ± 0.33 ± 0.25 
2.01 ± 0.13 ± 0.22 


0.28 ± 0.06 ± 0.05 % 
0.52 ± 0.18 ± 0.12 % 
0.85 ± 0.15 ± 0.19 % 


TABLE XL Average raw dilution for three flavor taggers in the electron and muon data 
corrected for wrong sign sequential decay trigger leptons and the bb to cc ratio. 


Flavor Tag 


(-^raw) 


Correction 




(z);,„(e))/(z);,„(/z)) 


JCSV (e) 
JCSV ifi) 


2.7 ± 0.5 % 
2.3 ± 0.4 % 


1.16 
1.26 


3.1 ± 0.6 % 
2.9 ± 0.5 % 


1.07 ± 0.28 


JCDV (e) 
JCDV (/i) 


12.1 ± 1.1 % 
9.7 ± 1.0 % 


1.15 
1.24 


14.2 ± 1.3 % 
12.0 ± 1.2 % 


1.18 ± 0.16 


SLT (e) 
SLT In) 


22.1 ± 1.3 % 

18.2 ± 1.3 % 


1.15 
1.24 


25.5 ± 1.5 % 
22.5 ± 1.6 % 


1.13 ± 0.10 
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FIG. 1. Jets distributions from the inclusive lepton data samples after the trigger lepton has 
been associated with a secondary vertex. The upper figure shows the number of jets with transverse 
momentum > 5 GeV/c per event. Approximately 60% of the events have more than one jet. 
The lower figure shows the separation in azimuth of the trigger-lepton jet and the jet with the 
largest pt in these events. 
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FIG. 2. The sample composition fits using p^^'- The left hand plots are for the electron data 
and the right hand plots are for the muon data. The fit values for the fraction of events from hh 
production are given in Table |. 
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FIG. 3. The sample composition fits using rricx- The left hand plots are for the electron data 
and the right hand plots are for the muon data. The fit values for the fraction of events from hh 
production are given in Table |. 
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FIG. 4. The observed momentum fraction, which is defined as VtIVt- The top plot shows two 
ranges of cluster mass (m'^^). The bottom plot shows two ranges of cluster transverse momentum 
(p^). The plots were obtained from the e-trigger 65 Monte Carlo sample. 
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FIG. 5. The reconstructed proper decay length ct for the e and trigger data. The pomts are 
the data. The soHd histogram is a combination of hb and cc Monte Carlo with the relative fractions 
given in Table |. All distributions are normalized to unit area. The upper plots show the data on 
a linear scale while the lower plots show it on a logarithmic scale. 
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FIG. 6. Distributions of soft-e and soft-/i p^"^ for e-trigger and ^-trigger events. The solid 
(dashed) distributions are for events where the charge of the soft and trigger leptons have opposite 
(hke) sign. The bin below is for events where the soft lepton was isolated. 
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FIG. 7. Soft-e and soft-/U tag raw dilution as a function of soft-lepton pq^^ for the e and 
trigger data. The sohd curve is the parameterization used for calculating the expected raw dilution 
on an event-by-event basis based on the soft-lepton p^?'. The dashed curves show the variation 
used in the evaluation of the systematic uncertainty on the raw dilution parameterization. 
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FIG. 8. Jet charge distributions for double (top) and single (bottom) vertex events for the e 
and /U trigger data. The data have been divided into events with a positively or negatively-charged 
trigger lepton (£+ or i^). Note that the separation between the ^+ and i'^ distributions is signifi- 
cantly larger for the double vertex events. 
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FIG. 9. The raw dilution Draw as a function of absolute jet charge \Qjet\ for double and single 
vertex events in the e and /x trigger data. Note the significantly higher dilution for double vertex 
events. 
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FIG. 10. The transverse decay length resolution functions, defined as the difference between 
the measured and true L^y in the bb Monte Carlo, for the e and /j, trigger data. The distributions 
were fit to three Gaussians and two exponentials. 
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FIG. 11. The same sign fraction as a function of the proper decay length ct. A representation 
of the unbinned maximum likehhood fit for Am^ and the Njj factors is superimposed on the data. 
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FIG. 12. Fit results for 400 fast Monte Carlo samples, simulating the e-trigger, SLT flavor 
tagged data. The top row of plots show distributions of the fitted values with arrows indicating the 
generated values for the samples. The middle plots show distributions of the statistical uncertainties 
on the fitted values with arrows indicating the statistical uncertainty for the fit to the real data. 
The bottom plots show distributions of the deviation from the true value, divided by the statistical 
uncertainty (pull distributions). The pull distributions have zero mean and unit width, as expected. 
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